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Abstract: In this work, we focus on the collapse of a spherically symmetric perturbation, with 
a classical top-hat profile, to study the nonlinear evolution of only viscous modified Chaplygin gas 
(VMCG) perturbations in Einstein’s gravity as well as in loop quantum Cosmology (LQC). In the 
perturbed region, we have investigated the natures of equation of state parameter, square speed of 
sound and another perturbed quantities. The results have been analyzed by numerical and graphical 
investigations. 


I. INTRODUCTION 


Recent years have witnessed the emergence of the idea 
of accelerating Universe and due to some observational 
results pj, 3 it is now established that the Universe is ac¬ 
celerating. This acceleration is caused by some unknown 
matter dubbed as dark energy. This dark energy has the 
positive energy density and strong negative pressure that 
satisfies equation of state parameter w = p/p < —1/3. 
The present acceleration is also confirmed by other ob¬ 
servations like large scale structure @, CMBR 0 and 
WMAP P, 0. The observations predict that the present 
Universe occupies ~ 4% ordinary matter, ~ 74% dark en¬ 
ergy and ~ 22% dark matter. The most simple candidate 
of dark energy is the cosmological constant A which sat¬ 
isfies the EoS parameter w = — 1 Q. Another candidates 
of dark energy are quintessence (where EoS parameter 
satisfies — 1 < w < —1/3) Q and phantom (where EoS 
parameter satisfies w < — 1) @. A brief review of dark 
energy models is found in the ref. 0- The unified dark 
fluid (UDF) model 0 0 was investigated extensively 
in the recent years. The main features of the UDF model 
are that it combines cold dark matter and dark energy 
and that it behaves like the cold dark matter and the dark 
energy at an early epoch and a late time respectively how¬ 
ever the approach remains purely phenomenological. In 
recent years, an increasing number of cosmological obser¬ 
vations suggest that our universe is filled with imperfect 
fluid which contains bulk viscosity which has negative 
pressure 0 , 0 . The viscous generalized Chaplygin gas 
model is one of the suitable candidate of unified dark en¬ 
ergy and cold dark matter as a unique imperfect dark 
fluid 0, 0 | . 

While dark energy may be modeled as a fluid with neg¬ 
ative pressure acting against gravitational collapse, dark 
matter (in its cold version) is a dust-like fluid with no 
pressure, therefore enhancing the collapse of matter per¬ 
turbations 00 - It is very important to investigate the 
evolutions of density perturbations within realistic cos¬ 
mological model (for recent reviews, see for instance 0 
and references therein). Fabris et al (20) studied the evo¬ 
lution of density perturbations in a universe dominated 
by the Chaplygin gas. Their model shows the required 
density contrast observed in large scale structure of the 


universe, however their approach remained Newtonian. 
Similar investigations were performed by Carturan and 
Finelli for the generalized Chaplygin gas 0] (see also 
0 |). Fernandes et al 0 ] have investigated the non¬ 
linear collapse of generalized Chaplygin gas in the frame 
of spherical top-hat collapse (STHC). Recently Carames 
0 investigated the spherical collapse model using vis¬ 
cous generalized Chaplygin gas. Li et al 0,0 have 
extended the above work by considering bulk viscosity in 
the general Chaplygin gas model. Besides the parame¬ 
ter a, they have also analyzed the effect of bulk viscos¬ 
ity on the structure formation of the variable generalized 
Chaplygin gas model which has a spherically symmetric 
perturbation. 

In our work, we focus on the collapse of a spherically 
symmetric perturbation, with a classical top-hat pro¬ 
file, to study the nonlinear evolution of viscous modified 
Chaplygin gas perturbations in Einstein’s gravity as well 
as in loop quantum Cosmology in sections II and III re¬ 
spectively. The conclusion is present in the last section 
IV. 


II. SPHERICAL TOP-HAT COLLAPSE MODEL 
OF VISCOUS MODIFIED CHAPLYGIN GAS IN 
EINSTEIN’S GRAVITY 


A. Basic Equations 


We consider a flat Friedmann Robertson-Walker 
(FRW) universe described by the following metric 

ds 2 = -dt 2 + a 2 (t) [dr 2 + r 2 d0 2 + r 2 sin 2 ddcp 2 ] . (1) 

Here a(t) is the scale factor of the universe. We assume 
that the spacetime is filled with only one component fluid 
having a bulk viscosity. The Einstein’s field equations are 
given by 


and 


H 2 


8nG 


3 


P , 


( 2 ) 


H = -4?rG(p + p), (3) 


where H(= |) is the Hubble parameter, p is the energy 
density and the effective pressure p can be expressed as 
follows 
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p = p d + n, 


(4) 
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which is the sum of the equilibrium pressure pde and the 
bulk pressure II = —£u7 7 , where u 7 is the four velocity of 
the fluid and £ is the coefficient of bulk viscosity and is a 
function of energy density. The first attempts at creating 
a viscosity theory of relativistic fluids were executed by 
Eckart J2{| and Landau and Lifshitz [27j who considered 
only a first-order deviation from equilibrium. The bulk 
viscous pressure II is represented by the Eckart’s expres¬ 
sion which is proportional to the Hubble parameter H 
with proportionality factor identified as the bulk viscos¬ 
ity coefficient £ which is defined by £ = £o p v , £o and v 
are constants. For simplicity, choosing v = n can be 
written as (HI 


n = -3 ZoHy/p. (5) 

The Chaplygin gas is generally characterized by the 
equation of state (pressure is inversely proportional to 
the energy density): p = —A/p. This equation of state 
has an interesting connection with the d-branes which 
are expressed via Nambu-Goto action [2fj. It also enjoys 
connections with the Newtonian hydrodynamical equa¬ 
tions. Further the Eddington-Born-Infeld model can be 
seen as an affine connection version for the Chaplygin gas 
approach j28j. The Chaplygin gas has been extensively 
studied within the unified dark energy-dark matter mod¬ 
els as well [29}. Therefore we choose the equation of state 
for modified Chaplygin gas as [30] 

p d = Ap -(6) 

P 

where A, B and a are constants which are constrained by 
the astrophysical data (for recent constraints, see a 
The modified Chaplygin gas satisfactorily accommodates 
an accelerating phase followed by a matter dominated 
phase of the universe. It is also consistent with the ob¬ 
servational studies dealing with the large scale structure 
& Hence combining the equations (4) to (6), we get 

P = A P- ^ - 3£oH^p, (7) 

which is called viscous modified Chaplygin gas (VMCG) 
tdtil . Hence by choosing the above equation of state, we 
extend previous works dealing with viscous linear equa¬ 
tion of state LB = 0) and viscous generalized Chaplygin 
gas ( A = 0) [H, 01 • It may also be termed as a ‘viscous 
unified dark fluid’. The conservation equation is given by 

p + 3H(p + p) =0 (8) 


Now using equations (7) and (8), we obtain the solution 


B C 

1 + A - V3 £ 0 + a 3 ( 1 + a )( 1 + A ~^ &) 

where C is a constant of integration. Now using the red- 
shift formula z = - — 1 (choosing a 0 = 1), equation (9) 
can be re-written in the form: 




p{z) = po \A a + (1 - 4,)(1 + z) 3 ( 1 +“)( 1 +^-'/3f°) | 

( 10 ) 

where po is the present value of the density and A s = 
(1+ A-Jl; o) c + B with 0 < A, < 1 and 1 + A > v^o- 
Hence the Hubble parameter is obtained as 


H(z ) — Hq {4l s + (1 — 7L s )x 


(1 + z ) 3 ( 1 +«)( 1 + A -v / 3&) l 1+1 


( 11 ) 


where S2o = 8 g ^ 2 ° is the present value of dimension den¬ 
sity parameter and Ho ~ 72 km s~ 1 Mpc~ 1 . The equation 
of state parameter is given by 


= -= A — 


B 

-,1+a 


3£o H 


P Vp ’ 

while the adiabatic sound speed reads 


c 2 _ _ A 
° s ~ dp ~ A 


aB £ 0 H 


t 1 + 0! 


( 12 ) 


(13) 


Note that the Chaplygin gas model is a dynamical dark 
energy model as well and dark energy perturbations play 
crucial role in dynamical dark energy models as well |34| . 


B. Spherical top-hat collapse model 

The spherical collapse (SC) which provides a way to 
glimpse into the nonlinear regime of perturbation theory 
was introduced firstly by Gunn and Gott :33). The SC 
describes the evolution of a spherically symmetric pertur¬ 
bation embedded in a homogeneous background, which 
can be static, expanding or collapsing. One assumes a 
spherical ‘top-hat’ profile for the perturbed region, i.e., 
a spherically symmetric perturbation in some region of 
space with constant density. Following the assumption 
of a top-hat profile, namely the density perturbation is 
uniform throughout the collapse, so the evolution of per¬ 
turbation is only time-dependent. 

The perturbed quantities p c and p c are related to their 
background counterparts by p c = p + Sp and p c = p + 
5p. Now the perturbed equation of state w c is given by 


p + Sp w + c/,5 
p + Sp 1 + S 


1 / B 3£ 0 H\ & 

1 + 5 \ p l + a y/p ) 1 + <5 


(14) 


where 5 = ^ and perturbed square of the sound speed is 
given by 


2 5p 

c ' = Tp 


Pc 
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A)Vp 

(15) 

In the spherical top-hat collapse model, the background 
evolution equations are still in the forms 3 S3, 23]: 


p = -3H(p+p), (16) 


a 

a 


4t tG 


(p + 3 p) 


(17) 


For the perturbed region, the basic equations which 
depend on local quantities can be written as 


p c = -3 h(p c + p c ) 


(18) 
























(19) 
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where h = - is the local expansion rate and r is the local 
scale factor and furthermore, h relates to local expansion 
rate in the STHC framework [23] . [35| , 

h = H + l ( 20 ) 

where 6 = S7 .v is the divergence of the peculiar velocity 

v. 

After some calculations [23|, [24[ , the dynamical equa¬ 
tions of density contrast S and 6 can be obtained in the 
following forms: 

S' = -~(c 2 e - w)S - [1 + w + (1 + c 2 e )S]-^— , (21) 

a cl ±1 


9 ’ 


9 9 2 

a 3 a 2 H 


—ns(i + 3c 2 e ) 


( 22 ) 


where, Q = S ^ 2 P and ' represents the derivative w.r.t. 
scale factor a. The above two equations can be re-written 
in the forms: 


d5_ 

dz 


3 {cl - w)S 
1 + z 


+ [1 + w + (1 + c 2 )S] 


H(z) ’ 


(23) 


d9 

dz 


9 9 2 

1 + z + 3 H{z) 


3 H{z) 
2(1 + z) 2 


nj(i + 3t^). 


(24) 


We have drawn the time varying parameters 
S , 9 , w, w c , c 2 , c 2 , h vs redshift z in figures 1-7 respectively 
for Einstein’s gravity. In all the figures we have taken the 
values of the parameters A = 0.9, B = 0.6, C = 0.2, £ 0 = 
0.2, a = 0.5, Hq = 72, flo = 1- From fig.l we observe that 
the parameter <5 increases from 0 to the value ~ 1.8 as the 
redshift z decreases from 5 to — 1. Similar nature hap¬ 
pens for the quantity 9 (increases from 0 to ~ 3.3) which 
is shown in fig.2. The equation of state parameter w vs 
2 is drawn in figure 3. We see that w decreases from 0.15 
to —1. So the VMCG generates initially normal fluid 
and after certain stage, it generates quintessence dark 
energy but phantom barrier does not happen, but mean¬ 
while the perturbed equation of state parameter w c does 
not so change and its value ~ 0.9 all over the time. The 
square speed of sound c 2 for our VMCG fluid system first 
increases from 0.61056 for 2 ~ 2 and ultimately increases 
to the value 0.61059 at z ~ —0.5 (fig.5). But, perturbed 
square speed of sound c 2 actually increases from 0.74 (at 
2 ~ 5) to 0.9 after z ~ 4.8 (fig.6). Also the expansion 
rate h vs z has been drawn in figure 7. We observe that 
h first increases from 0 (from z ~ 5) to ~ 60.1 around 
z ~ 0.4 and after that suddenly decreases to zero (nearly 
* = -!)• 


III. SPHERICAL TOP-HAT COLLAPSE MODEL 
OF VISCOUS MODIFIED CHAPLYGIN GAS IN 
LOOP QUANTUM COSMOLOGY 
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Loop Quantum Gravity (LQG) is a canonical quantiza¬ 
tion of gravity based upon Ashtekar connection variables. 
LQG is an important frontier to explore the quantum 


Figs. 1-4: Plots of S,6,w,w c vs redshift z in Einstein’s 
gravity. 
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Fig-7 

Figs. 5-7: Plots of c 2 , c 2 , h vs redshift z in Einstein’s gravity. 


gravity effects in cosmology. Some of its implications in¬ 
cludes the prediction of cosmic inflation in the early uni¬ 
verse iil , late time cosmic acceleration [37] and primor¬ 
dial gravitational waves fUj. The field equations of LQC 
admit attractor solutions which are of enormous cosmo¬ 
logical interest. The cosmological perturbation theory 
within LQC has also been explored in [39]. We here 
consider the flat homogeneous and isotropic universe de¬ 
scribed by FRW metric, so the modified Einstein’s field 
equations in LQC are given by |"iTil llfjt 


H 2 


8nGp 



(25) 


and 


H = -A7rG(p + p) (l-^) , (26) 

where p\ = le7T 2^3G 2 h is called the critical loop quantum 
density, 7 is the dimensionless Barbero-Immirzi parame¬ 
ter. Now the fluid is considered as viscous modified Chap¬ 
lygin gas and in this case, the density is given in equation 
(10) and in LQC, the Hubble parameter is obtained as 


H(z ) = H 0 



x 



3 H 2 n 0 
Pi 



A s )( 1 + 2 )3(1+«)(1+A-V3 £o)| T F 


A s )(l + ^) 3 ( 1 +«)( 1 +^-v / 3«o)| 1+1 


(27) 

Similar to the above section, the equation (23) will be the 
same and equation (24) modifies to the form 

dff _ _e_ e 2 3 h{z) 

dz 1 + z 3 H(z) 2(1 + z) 2 


(1 + 3Cg) SnGpi ^ + + ^ c e) + ^ w ~ ^ 

( 28 ) 

In this case also, the equation of state w c and square of 
the sound speed are given in (14) and (15). We have 
drawn the time varying parameters 3, 9, w, w c , c 2 , c 2 ,h vs 
redshift z in figures 8-14 respectively for LQC. In all the 
figures we have taken the values of the parameters A = 
0.9, B = 0.6, C = 0.2, £ 0 = 0.2, a = 0.5, H 0 = 72, = 1. 

From fig.8 we observe that the parameter S decreases from 
a certain value (around z ~ 4). Also the quantity 9 (in¬ 
creases from 0 to ~ 4.4) during z = 5 to z = 4.8 and after 
that 9 decreases to zero (upto z ~ — 1) which is shown in 
fig.9. The equation of state parameter w vs z is drawn in 
figure 10 which shows that w decreases from some neg¬ 
ative value to —0.84 as z decreases. So the VMCG in 
our considered LQC model generates like quintessence 
dark energy but phantom barrier does not happen. In 
the meanwhile the perturbed equation of state parame¬ 
ter w c changes from some positive value (< 1) to ~ —0.9 
(fig.11). The square speed of sound c 2 for our VMCG 
fluid system first decreases to the value 0.61 at z ~ —0.5 
(fig.12). But, perturbed square speed of sound c 2 actu¬ 
ally increases from some fractional value (< 0.5) to ~ 0.02 
(fig.13). Also the expansion rate h vs z has been drawn 
in figure 14. We observe that h first increases from some 
positive value (from z ~ 5) to ~ 4 around z ~ 4.8 and 
after that it decreases to zero (nearly z = —1). 


IV. CONCLUSIONS 

In this work, we mainly focused on the collapse of 
a spherically symmetric perturbation, with a classical 
top-hat profile, to study the nonlinear evolution of only 
viscous modified Chaplygin gas (VMCG) perturbations 
in Einstein’s gravity as well as in loop quantum Cosmol¬ 
ogy (LQC). The background model is considered as flat 
FRW metric. Since we know that modified Chaplygin 
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Figs. 8-11: Plots of 5,0,w,w c vs redshift z in LQC. 
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Figs. 12-14: Plots of c^, c^, h vs redshift z in LQC. 


gas (MCG) is the unified model of dark matter and 
dark energy. So we have not assumed any dark matter 
separately. In this occasion, we have assumed viscous 
modified Chaplygin gas (VMCG) (which is also the 
unified model as established) by including the viscosity 
term in the equation of state in MCG. The density and 
the Hubble parameter have been calculated in terms 
of redshift z in both gravity models. The equation 
of state parameter, square speed of sound have also 
been calculated for our considered VMCG model in 
Einstein’s gravity and LQC also. Next step, the density 
perturbation for our top-hat collapsing profile has 
been investigated. In the perturbed region, we have 
investigated the natures of equation of state parameter, 
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square speed of sound and another perturbed quantities 
like 5, 9 , h etc. The dynamical equations of density 
contrast S and 6 have been found in both gravity models. 
Analytically, it is not possible to find out the natures of 
the perturbed quantities. So numerically and graphically 
we analyzed the natures of the perturbed quantities 
which are given as follows: 

• Einstein's Gravity. We have drawn the time 
varying parameters S,9,w,w c ,c2 : c^,,h vs redshift z in 
figures 1-7 respectively for Einstein’s gravity. In all the 
figures we have taken the values of the parameters A = 
0.9, B = 0.6, C = 0.2, £ 0 = 0.2, a = 0.5, H 0 = 72, Q 0 = 1. 
From hg.l we observe that the parameter <5 increases 
from 0 to the value ~ 1.8 as the redshift 2 decreases 
from 5 to — 1. Similar nature happens for the quantity 6 
(increases from 0 to ~ 3.3) which is shown in fig.2. The 
equation of state parameter w vs z is drawn in figure 3. 
We see that w decreases from 0.15 to —1. So the VMCG 
generates initially normal fluid and after certain stage, it 
generates quintessence dark energy but phantom barrier 
does not happen, but meanwhile the perturbed equation 
of state parameter w c does not so change and its value 
~ 0.9 all over the time. The square speed of sound (? s for 
our VMCG fluid system first increases from 0.61056 for 
2 ~ 2 and ultimately increases to the value 0.61059 at 
2 ~ —0.5 (fig.5). But, perturbed square speed of sound 
Cg actually increases from 0.74 (at z ~ 5) to 0.9 after 
2 ~ 4.8 (fig.6). Also the expansion rate h vs z has been 
drawn in figure 7. We observe that h first increases from 
0 (from z ~ 5) to ~ 60.1 around z ~ 0.4 and after that 
suddenly decreases to zero (nearly z = —1). 


• LQC: We have drawn the time varying parameters 
6,9,w,w c ,c%,c%,h vs redshift z in figures 8-14 respec¬ 
tively for LQC. In all the figures we have taken the 
values of the parameters A = 0.9, B = 0.6, C = 0.2, £ 0 = 
0.2, a = 0.5, Hq = 72, Q 0 = 1- From fig.8 we observe that 
the parameter S decreases from a certain value (around 
z ~ 4). Also the quantity 9 (increases from 0 to ~ 4.4) 
during z = 5to2: = 4.8 and after that 8 decreases to zero 
(upto 2 ~ — 1) which is shown in fig.9. The equation 
of state parameter w vs z is drawn in figure 10 which 
shows that w decreases from some negative value to 
—0.84 as 2 decreases. So the VMCG in our considered 
LQC model generates like quintessence dark energy but 
phantom barrier does not happen. In the meanwhile 
the perturbed equation of state parameter w c changes 
from some positive value (< 1) to ~ —0.9 (fig.11). The 
square speed of sound (? s for our VMCG fluid system first 
decreases to the value 0.61 at 2 ~ —0.5 (fig.12). But, 
perturbed square speed of sound Cg actually increases 
from some fractional value (< 0.5) to ~ 0.02 (fig.13). 
Also the expansion rate ft vs 2 has been drawn in figure 
14. We observe that h first increases from some positive 
value (from 2 ~ 5) to ~ 4 around 2 ~ 4.8 and after that 
it decreases to zero (nearly 2 = — 1). 
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